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NOMENCLATURE 


“ef 
b ,  c 

{F} 

[F] 

F,, 

f(t) 

[I] 

i 

[M] 

Q’ 


in 


It 


to  ,  (0, 
a  b 

to 

n 

Qk 

{q(t)} 

{q(t)} 

{q(t)} 

[4>] 

[♦] 

‘•’e.k 

{♦}e 

{x(t)} 


x’j^Cto) 


orthogonality  between  modes  e  and  f 
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Subscripts 

e  =  mode  number 

f  =  mode  number 

i  =  iteration  number 

k  =  mode  number 

8.  =  physical  coordinate  number 

m  =  counter  for  physical  coordinates 

r  =  number  of  physical  coordinates 

t  =  target  mode,  mode  to  be  measured 
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I.  INTRODUCTION 


Formulation  of  accurate  dynamic  models  of  complex  structures  requires  the 
experimental  measurement  of  the  structure's  natural  modes  of  vibration.  A 
large  number  of  test  procedures  and  methodologies  have  been  developed  over  the 
years  to  reduce  the  complexity  and  cost  of  establishing  this  data.  These 
procedures  can  be  categorized  into  those  that  measure  the  dynamic  properties 
directly  (e.g..  Refs.  1-5)  and  those  that  identify  the  needed  properties  from 
curve  fitting  of  frequency  response  functions  or  time  domain  free  response 
data  (e.g..  Refs.  6-15).  The  decision  as  to  which  procedure  to  use  is  highly 
dependent  on  the  test  conditions  and  requirements. 

The  direct  measurement  procedures  have  several  advantages.  However,  the 
principal  advantage  is  that  the  dynamic  properties  are  measured  directly. 

Thus,  the  disadvantage  of  interpreting  numerical  curve  fitting  results,  as 
required  by  the  analytical  mode  identification  procedures,  is  avoided.  The 
principal  disadvantage  of  the  direct  measurement  procedures  is  that  they  are 
relatively  time  consuming. 

The  most  widely  used  direct  measurement  mode  survey  test  procedure  is  the 
multi-shaker  sine-dwell  approach.  The  procedure  consists  of  using  one  or  more 
shakers  to  exert  sinusoidally  varying  forces  on  the  test  article.  The 
frequency  of  excitation  and  the  shaker  force  levels  are  adjusted  to  best  iso¬ 
late  the  target  mode  response  from  all  other  modes.  Mode  parameters  are  then 
established  from  direct  measurement  of  the  forced  vibration. 

Available  sine-dwell  test  procedures  that  establish  mode  parameters  by 
direct  measurement  require  iterative  adjustment  of  the  shaker  force  levels  and 
excitation  frequency.  The  adjustment  process  is  particularly  time  consuming 
because  shaker/structure  interaction  creates  a  need  to  adjust  phase  between 
the  shaker  force  control  signals  to  maintain  colinear  forces.  This  is  further 
complicated  when  performing  multi-shaker  sine-sweeps  since  the  phase  between 
shaker  forces  will  vary  as  a  function  of  test  article  response.  Frequency 
sweep  results  obtained  without  maintaining  proper  phasing  of  the  forces  are 
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usually  misleading  and  not  indicative  of  the  degree  of  mode  isolation.  An 
alternative  is  to  perform  a  series  of  sine  dwells  and  adjust  the  phase  at  each 
frequency.  Obviously,  this  is  very  time  consuming. 

It  is  the  purpose  of  this  discussion  to  introduce  a  more  efficient 
approach  to  multi-shaker  sine-dwell  mode  survey  testing.  The  procedure  derives 
from  the  principle  that  the  total  response  of  a  linear,  elastic  structure  con¬ 
sists  of  the  superposition  of  responses  to  the  individual  excitation  sources. 
Assuming  that  linearity  has  been  confirmed  experimentally,  frequency  response 
functions  obtained  with  single  shakers  can  be  scaled  and  summed  numerically  to 
yield  frequency  response  functions  corresponding  to  multiple-shaker  excitation. 
Thus,  the  adjustments  to  establish  force  levels  needed  to  isolate  each  mode 
for  measurement  can  be  performed  numerically.  The  final  forces  can  then  be 
applied  to  the  structure  and  the  entire  mode  vector  measured. 

This  also  introduces  the  possibility  of  identifying  mode  parameters  from 
single-shaker  frequency  response  functions  using  traditional  multi-shaker 
sine-dwell  testing  techniques.  To  accomplish  this,  all  that  is  required  is 
that  frequency  response  functions,  corresponding  to  each  excitation  location, 
be  measured  for  all  accelerometer  locations.  Once  the  appropriate  force 
levels  for  mode  isolation  have  been  established,  all  corresponding  frequency 
response  functions  can  be  scaled  and  summed  numerically.  The  mode  parameters 
can  then  be  obtained  directly  from  the  resulting,  multi-shaker  equivalent, 
frequency  response  functions. 
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II.  RESPONSE  TO  HARMONIC  EXCITATION 


The  behavior  of  a  large  class  of  linear,  elastic  structures  subjected  to 
multiple  harmonic  forces  {F}f(t)  can  be  described  in  modal  coordinates  by 
the  matrix  differential  equation  of  motion 

[I]{q(t)}  +  [2Cm^]{q(t)}  +  [tO^]{q(t)}  =  [<j)  ]  ^{F}  f  ( t  )  (1) 

The  coordinate  transformation  between  physical  coordinates  {x(t)}  and  modal 
coordinates  {q(t)} 

{x(t)}  =  [(|)]{q(t)}  (2) 

and  the  matrix  of  mode  shape  vectors  [(^]  has  been  normalized  with  respect  to 
the  mass  matrix  such  that 


[4>]'^[M]t4>]  =  [I] 


(3) 


Equation  (1)  consists  of  uncoupled  second  order  differential  equations  of 
the  form 


■4<t) 


*  "nk'lk 


(t)  =  J  4. 

U1 


(4) 


By  solving  Equation  (4)  for  each  mode,  the  frequency  response  function  for  a 
typical  physical  coordinate  can  be  obtained: 


where 
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In  Equation  (5)  the  real  part  corresponds  to  the  component  of  acceleration 
response  that  is  colinear  with  the  reference  force  time  history,  i.e.,  the 
coincident  (Co)  response.  The  imaginary  part  corresponds  to  the  component  of 
acceleration  response  that  is  90®  out  of  phase  with  the  reference  force  time 
history,  i.e.,  the  quadrature  (Qd)  response.  The  term  in  Equation  (5)  is 
the  kth  mode  modal  force  amplitude.  The  solution  to  Equation  (4)  was  obtained 
by  establishing  the  particular  solution  corresponding  to  the  modal  force 
Equation  (4),  however,  can  also  be  solved  by  superposition  of 
particular  solutions  corresponding  to  the  individual  terms  that  comprise 
Using  this  approach,  the  frequency  response  function  for  a 
typical  physical  coordinate  is 


Xg^(co)/f  (w) 


I  {  I 

m=l  (k=l 


+  iQ, 


><)>  I 
mk  m  1 


(10) 


Equations  (5)  and  (10)  yield  equivalent  frequency  response  functions.  How¬ 
ever,  by  using  Equation  (10)  the  frequency  response  functions  corresponding  to 
each  exciter  are  established  individually.  The  superposition  of  these  func¬ 
tions  then  yields  the  frequency  response  function  corresponding  to  multiple 
shaker  excitation.  The  advantage  of  using  Equation  (10)  is  that  the  response 
to  a  different  set  of  force  levels  can  be  obtained  by  simply  scaling  the 
appropriate  transfer  functions  and  then  summing  them. 


This  suggests  an  efficient  approach  to  multi-shaker  sine-dwell  mode  survey¬ 
ing  and  introduces  the  possibility  of  using  these  techniques  in  a  post-test 
mode  identification  procedure.  Single-shaker  excitation  could  be  used  to 
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obtain  frequency  response  functions  for  various  locations  of  interest.  The 
shaker  can  then  be  moved  to  another  location  (or  another  shaker  can  be  used) 
and  a  new  set  of  frequency  response  functions  obtained.  This  can  be  repeated 
for  several  shaker  locations.  Now,  to  obtain  the  frequency  response  functions 
corresponding  to  multiple  shaker  excitation,  all  that  is  needed  is  that  the 
appropriate  functions  be  summed.  Also,  if  different  force  levels  are  desired, 
the  corresponding  frequency  response  functions  need  only  to  be  scaled  before 
being  added. 


11 


12 


III.  ISOLATION  OF  TARGET  MODE 


Before  introducing  the  mode  identification  procedure,  a  systematic 
approach  to  establishing  the  required  force  levels  for  mode  isolation  will  be 
described.  However,  practical  aspects  of  mode  shape  modal  contamination  need 
to  be  addressed  first. 

Assume  for  the  moment  that  the  test  article  can  be  harmonically  excited  at 
all  its  degrees  of  freedom.  Equation  (5)  indicates  that  the  test  article  can 
be  made  to  vibrate  in  a  single  normal  mode  by  appropriately  adjusting  the 
shaker  force  levels  such  that  the  only  nonzero  modal  force  is  that  of 
the  target  mode.  In  practice,  however,  the  number  of  degrees  of  freedom 
greatly  exceeds  the  number  of  available  shakers.  Therefore,  perfect  isolation 
is  not  possible  and  some  degree  of  modal  contamination  will  have  to  be 
accepted. 

Modal  contamination  can  be  quantified  by  calculating  the  mass  weighted 
orthogonality  between  measured  mode  vectors  that  have  been  normalized  to  unit 
modal  mass.  Let  the  "measured"  mode  vectors  and  consist  of 

the  quadrature  components  of  total  response  and  be  defined  as 

+  b{4>}^  (11) 

and 

{r}f  =  {4>}f  +  c{(j>}^  (12) 

The  orthogonality  between  and  {<!>"'} ^  is  therefore 
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(14) 


/U  2  .2  ^,2  2,1/2 

=  (b  +  c)/(l  +  c  +  b  +  b  c  ) 

=  b  +  c  for  small  b  and  c  (15) 

It  should  be  noted  that  contamination  by  additional  modes  would  only  introduce 
higher  order  terms  in  Equation  (14).  Thus,  to  a  first  approximation.  Equation 
(15)  would  still  provide  the  mass  weighted  orthogonality  between  and  • 

It  is  generally  accepted  in  industry  that  mode  shapes  of  acceptable 
quality  have  been  measured  if  is  less  than  0.10.  This  requirement 
will  be  satisfied  if 


(Q^/  )(Qf/Qg)  <  b  (16) 

(fi^V  Qj  )(Qg/Qf)  <  c  (17) 

and 

b  +  c  <  0.10  (18) 

The  above  equations  indicate  that  adequate  mode  isolation,  from  an  off- 

resonance  mode,  will  exist  if  the  product  of  the  quadrature  admittance  ratio 

and  modal  force  ratio  is  less  than  b  when  {(b"’}  is  measured  and  less  than  c 
m  e 

when  is  measured. 


For  lightly  damped  structures  the  ratio  (nj^/Q^.)  will  decrease  rapidly 
with  increasing  frequency  separation  between  the  two  modes.  Thus,  careful 
attention  to  the  applied  forces,  which  determine  the  modal  force  ratio  (Q]j/Q^)i 
is  necessary  only  for  modes  close  in  frequency  to  the  target  mode.  Since  this 
is  usually  a  small  number,  only  a  few  shakers  will  be  needed  to  satisfy  the 
requirements  of  Equations  (16),  (17),  and  (18).  Therefore,  we  shall  proceed 
with  the  understanding  that  perfect  isolation  of  a  target  mode  is  not  neces¬ 
sary,  and  that  adequate  isolation  can  be  achieved  with  a  relatively  small 
number  of  shakers. 
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In  Reference  4,  Anderson  formulated  a  systematic  approach  to  isolate  and 
measure  modes  using  multi-shaker  sine-dwell  excitation.  Anderson's  procedure 
recognizes  the  natural  selectivity  of  lightly  damped  jtructui.  i  to  greatly 
decrease  their  quadrature  response  with  increasing  separation  between  the 
natural  frequency  and  excitation  frequency.  Thus,  as  indicated  by  Equations 
(16)  and  (17),  careful  attention  to  the  applied  forces  will  generally  be 
necessary  only  for  modes  in  close  frequency  proximity  to  the  target  mode. 

Data  presented  in  Reference  4  and  practical  test  experience  indicate  that 
for  lightly  damped  structures,  generally  only  modes  with  frequencies  within 
10  to  15%  of  the  target  mode  frequency  need  to  be  considered  when  measuring 
the  target  mode.  To  isolate  the  target  mode  response  from  these  contaminating 
modes,  Anderson  proposed  the  selective  orthogonal  excitation  (SOREX)  test 
logic.  The  procedure  consists  of  using  successive  estimates  of  mode  shape 
values  to  establish  force  vectors  whose  corresponding  target  mode  modal  force 
is  much  larger  than  that  of  the  modes  close  to  it  in  frequency. 

This  is  accomplished  by  establishing  isolation  groups  of  one  or  more  target 
modes  and  all  other  modes  within  10  to  15%  in  frequency.  Shaker  locations 
equal  to  the  number  of  modes  in  the  isolation  group  are  then  selected.  Using 
single  shaker  excitation  frequency  response  functions  at  the  shaker  locations 
are  measured.  From  these  transfer  functions,  estimates  of  shaker  location 
mode  shape  values  are  obtained.  These  approximate  mode  shape  amplitudes  are 
then  used  to  calculate  the  relative  force  levels  to  be  applied  with  multiple 
shakers ,  i .e . , 


[F](i)  =  ^  (19) 

The  kth  column  of  [F]q)  consists  of  the  relative  force  levels  to  be  applied 
by  each  shaker  to  improve  the  isolation  of  the  kth  mode.  The  superscript  (i) 
designates  the  iteration  number.  The  columns  of  consist  of  the  esti¬ 

mated  mode  shape  values  (quadrature  components)  for  each  shaker  location.  The 
kth  column  of  corresponds  to  the  kth  mode  in  the  isolation  group. 
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Each  calculated  force  vector  is  then  applied  to  the  structure,  the  fre¬ 
quency  of  the  excitation  is  adjusted  to  maximize  quadrature  response,  and  new 
mode  shape  values  at  the  shaker  locations  are  measured.  These  values  are  then 
used  in  Equation  (19)  to  calculate  a  set  of  more  refined  forces.  Typically, 
only  one  or  two  iterations  are  needed  to  obtain  adequate  isolation.  Once  this 
is  achieved,  the  entire  mode  shape  is  measured. 

Should  an  isolation  group  consist  of  more  modes  than  there  are  available 
shakers,  force  patterns  can  still  be  obtained  using  a  least-square  error 
approximation: 

[Fl(i)  . 

However,  the  isolation  of  target  modes  will  generally  not  be  as  complete  as 
when  the  number  of  shakers  equals  the  number  of  modes  in  the  isolation  group. 
An  example  is  presented  in  Reference  A. 
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IV.  MODE  IDENTIFICATION  PROCEDURE 


A  more  efficient  approach  to  multi-shaker  sine-dwell  testing  than  prac¬ 
ticed  heretofore  will  now  be  proposed.  We  will  take  the  perspective  that  the 
procedure  will  be  implemented  on  a  small  laboratory  computer.  It  will  be 
assumed  that  standard  data  acquisition  and  Fast  Fourier  Transform  software  are 
available  as  part  of  the  total  mode  identification  package  to  be  described. 
Furthermore,  we  will  assume  that  sufficient  accelerometers  have  been  deployed 
such  that  motions  of  all  resonant  components,  of  all  modes  within  the  frequency 
range  of  interest,  can  be  measured. 

The  first  step  of  the  procedure  is  to  acquire  the  needed  frequency 
response  functions.  We  begin  by  selecting  a  number  of  potential  shaker  loca¬ 
tions.  Then,  using  random  or  slow  sinusoidal  sweep  excitation,  we  obtain 
frequency  response  functions  associated  with  single-shaker  forcing  at  each  of 
these  locations.  It  should  be  noted  that  it  is  also  possible  to  obtain 
simultaneously  several  sets  of  frequency  response  functions,  each  associated 
with  single-shaker  excitation,  by  forcing  the  test  article  with  multi-shaker 
uncorrelated  random  excitation  (Refs.  16  and  17).  If  post-test  mode 
identification  is  desired,  frequency  response  functions  for  all  accelerometers 
must  be  acquired.  Otherwise,  at  this  time,  data  are  needed  only  for  the 
potential  multiple  shaker  locations  and  any  other  coordinates  whose  response 
will  aid  in  identifying  all  modes  in  the  frequency  range  of  interest. 

These  frequency  response  functions  are  then  reviewed  to  determine  pre¬ 
liminary  natural  frequencies.  Mode  isolation  groups  are  established  consisting 
of  one  or  more  target  modes  and  any  other  modes  sufficiently  close  in  frequency 
to  warrant  inclusion  in  determining  multiple  excitation  force  levels.  As  dis¬ 
cussed  previously  all  modes  within  10  to  15%  in  frequency  to  the  target  modes 
should  generally  be  included. 

From  the  available  shaker  locations,  those  locations  with  strong  resonances 
and  phase  reversals  are  selected  for  each  isolation  group.  The  number  of 
shaker  locations  selected  should  equal  the  number  of  modes  in  the  isolation 
group.  Restricting  our  attention  to  one  isolation  group,  we  obtain  for  each 
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mode  in  the  group  estimates  of  mode  shape  values  at  the  selected  shaker  loca¬ 
tions  and  form  Then,  we  use  Equation  (19)  to  calculate  the  first 

iteration  force  levels  If  fewer  shaker  locations  than  modes  in  an 

isolation  group  are  available.  Equation  (20)  must  be  used. 

The  columns  of  are  the  force  levels  needed  to  better  isolate  each 

mode  in  the  isolation  group.  Now,  rather  than  applying  these  multiple  force 
levels  to  the  test  article  and  measuring  new  estimates  of  shaker  location  mode 
shape  values,  we  will  take  advantage  of  Equation  (10).  For  each  mode  in  the 
isolation  group,  the  following  steps  must  be  followed.  First,  we  use  the  force 
levels  from  the  appropriate  column  of  [F]^^  to  scale  the  corresponding 
transfer  functions.  Then,  for  each  coordinate  we  add  these  scaled  frequency 
response  functions. 

For  example,  assume  two  shakers  at  locations  0.  and  fi.+l  were  selected,  and 
the  first  iteration  force  levels,  for  the  first  mode  in  the  isolation  group, 
are  F^^  and  p  respectively.  Then,  all  frequency  response  functions 

corresponding  to  single-shaker  excitation  at  location  ft  must  be  multiplied 
by  Fg^  Likewise,  all  frequency  response  functions  associated  with 
single-shaker  excitation  at  location  ft+1  must  be  scaled  by  2^*  Finally, 

for  each  coordinate  the  scaled  frequency  response  functions  are  added.  The 
resulting  functions  are  equivalent  to  those  that  would  be  measured  if  the  test 
article  were  excited  simultaneously  with  and  2^* 

These  multi-shaker  equivalent  frequency  response  functions  can  then  be 
reviewed  to  determine  if  adequate  isolation  has  been  achieved.  If  isolation 
is  not  adequate,  better  estimates  of  shaker  location  mode  shape  values  can  be 
obtained  from  the  new  transfer  functions.  The  above  steps  are  repeated  for 
each  mode  in  the  isolation  group.  This  will  yield  new  estimates  for 
subsequently  a  new  set  of  more  refined  forces  [Fj^^. 
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Typically,  only  one  or  two  iterations  will  be  required  to  obtain  the 

forces  that  will  produce  adequate  isolation.  The  refined  forces  can  then  be 

applied  to  the  structure,  the  frequency  of  excitation  adjusted  to  maximize 

target  mode  quadrature  response,  and  the  entire  mode  vector  measured.  If 

frequency  response  functions  were  obtained  for  all  accelerometer  locations, 

the  above  scaling  and  summing  operations  could  be  performed  for  all  coordinate 

locations.  The  mode  shape  is  then  obtained  as  the  peak  quadrature  response, 

at  0)  ,  from  each  transfer  function.  Thus,  each  mode  shape  can  be  isolated 

nt 

and  recorded  without  collecting  any  additional  data  off  the  structure. 

Damping  should  be  obtained  from  the  coincident  response  functions  estab¬ 
lished  in  the  final  iteration  using 

C  =  ((/  -  J)/h  co^  (21 

t  a  b  nt 

In  the  above  equation  is  the  critical  damping  ratio,  is  the  target  mode 

natural  frequency,  and  and  (co^  ^  “b^  frequencies  adjacent  to 

at  which  the  coincident  response  peaks. 
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V.  NUMERICAL  SIMULATION  OF  A  TEST 


To  demonstrate  the  mode  identification  procedure,  an  analytical  simulation 
of  a  test  will  be  used.  The  "test  article"  mass  and  stiffness  properties  are 
presented  in  Table  1.  The  frequency  response  functions  for  single-shaker  exci¬ 
tation  were  obtained  using  the  closed  form  complex  solution  to  the  equations 
of  motion  [Eq.  (5)].  Critical  damping  ratios  of  0.02  were  assigned  to  each  of 
the  eight  modes.  Response  values  were  calculated  every  0.60  rad/sec. 

Table  1.  Test  Article  Mass  and  Upper  Triangular  Stiffness 
Matrix  Elements 


i.  j 

k. 

i,J 

M. 

i,  j 

k. 

i,J 

M. 

1,1 

18 

0.00055 

5,5 

6100 

0.055 

1.2 

-18 

5,7 

-200 

CM 

CM 

4118 

0.055 

6,6 

4018 

0.055 

2,3 

-  4000 

6,7 

-  4000 

2,4 

-100 

6,8 

-18 

3,3 

6100 

0.055 

7,7 

6202 

0.055 

3,5 

-100 

7,8 

-2 

4,4 

4100 

0.055 

8,8 

20 

0.00055 

4,5 

-  4000 

IK/- 

Ib/in. 

(b)  ,,  2,. 

lb-sec  /in. 

We  begin  the  simulated  test  by  exciting  the  structure  with  a  single  shaker 
at  coordinate  and  then  at  coordinate  X^.  We  process  the  response  data 
at  various  coordinates  into  frequency  response  functions.  Two  of  these  trans¬ 
fer  functions  are  presented  in  Figure  1.  A  review  of  this  data  indicates  the 
presence  of  three  modes  in  the  19  to  23  Hz  range  and  two  modes  in  the  28  to 
32  Hz  range.  In  this  simulation  we  will  restrict  our  attention  to  modes  with 
frequencies  below  50  Hz. 
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D0F1  CO/QD  RESP 


Figure  1.  Single-Shaker  Excitation  Frequency  Response  Functions, 
Coincident  Component  -  -  Quadrature  Component  _ 

We  will  need  two  isolation  groups.  The  first  group  will  consist  of  the 
three  modes  in  the  19  to  23  Hz  range,  and  the  second  group  will  consist  of  the 
two  modes  in  the  28  to  32  Hz  range.  Since  the  first  group  consists  of  three 
modes,  three  excitation  locations  must  be  selected.  Data  from  two  of  these 
locations  can  then  be  used  to  isolate  the  two  modes  in  the  second  group. 

We  select  coordinates  X.,  X  ,  and  X,  for  the  multiple  excitation  locations. 
We  start  by  forcing  the  structure  with  a  single  shaker  at  coordinate  X^,  and 
then  process  the  response  data  at  all  eight  coordinates  into  frequency  response 
functions.  The  transfer  functions  for  coordinates  X„,  X^.  and  X,  are  presented 
in  column  one  of  Figure  2.  We  then  excite  the  test  article  at  coordinate  X^ 
and  again  establish  frequency  response  functions  for  all  eight  coordinates. 

The  functions  for  coordinates  X  ,  X^  and  X,  are  presented  in  the  second 
column  of  Figure  2.  This  is  repeated  a  third  time  with  excitation  at  coordi¬ 
nate  X^.  The  corresponding  frequency  response  functions  for  coordinates 
X^,  X^  and  X^  are  presented  in  the  third  column  of  Figure  2. 
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Figure  2.  Single-Shaker  Coincident  (-  -  -)  and  Quadrature  ( - ) 

Frequency  Response  Functions  for  Degrees  of  Freedom  2,  5, 
and  6 


To  establish  the  required  multiple-shaker  force  levels,  estimates  of  mode 
shape  values  (quadrature  components)  at  the  shaker  locations  are  needed.  For 
the  first  isolation  group,  we  will  use  the  frequency  response  functions 
obtained  with  single-shaker  excitation  at  coordinate  .  The  approximate 
mode  shape  values  are  simply  the  peak  quadrature  response  near  the  estimated 
natural  frequency  of  each  mode.  If  an  obvious  peak  is  not  discernible,  the 
value  at  the  estimated  frequency  of  the  mode  can  be  u^ed. 


The  estimated  shaker  location  mode  shape  values  are  now  used  to  form  the 
truncated  mode  matrix  (The  subscript  denotes  the  iteration  number.) 

For  our  three-mode  isolation  group,  we  obtain 
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63.658 

71.351 

85.059 


-121.205 
-  10.375 
183.381 


49.642 

-57.977 

53.993 


By  substituting  the  above  matrix  into  Equation  (19),  the  multiple  excitation 
force  levels  needed  to  improve  the  isolation  of  each  mode  are  obtained. 

The  calculated  force  patterns,  scaled  to  a  maximum  value  of  unity,  are: 

0.644  -  1.0  0.635 

1.0  -  0.089  -  1.0 

0.482  0.823  0.363 

Next,  we  will  establish  the  frequency  response  functions  associated  with 
the  calculated,  multiple-shaker  force  levels.  We  begin  with  the  force  levels 
calculated  for  the  first  mode,  i.e.,  the  first  column  of  [F]^.  First,  we 
multiply  the  shaker  location  frequency  response  functions  obtained  with  single¬ 
shaker  excitation  at  coordinate  (first  column  of  Fig.  2)  ^2  1  “  0*6^^ 

(the  subscripts  2  and  1  denote  the  coordinate  and  mode  number,  respectively). 
Next,  we  multiply  the  functions  obtained  with  single-shaker  excitation  at 

coordinate  X  (second  column  of  Fig.  2)  and  coordinate  X  (third  column  of 
J  o 

Fig.  2)  by  Fc  ,  =  1.0  and  F  .  =  0.482,  respectively.  Finally,  for  each 
coordinate  (degree  of  freedom),  the  three  scaled  frequency  response  functions 
are  summed.  We  repeat  this  procedure  using  the  force  levels  calculated  for 
the  second  mode,  and  then  the  third  mode. 

The  resulting  frequency  response  functions  for  coordinate  X^  are  presented 
in  Figure  3  (frequency  response  functions  (b),  (c),  and  (d)).  These  response 
functions  are  those  that  would  be  obtained  if  the  multiple  force  levels  in  each 
of  the  three  columns  of  [F]j^  were  applied  to  the  test  article  one  set  at  a 
time  As  can  be  ascertained  from  frequency  response  function  (b),  isolation 
of  the  first  mode  relative  to  the  other  two  modes  has  improved  substantially. 
Likewise,  response  functions  (c)  and  (d)  indicate  that  the  second  and  third 
modes  are  also  much  better  isolated  by  the  corresponding  calculated  force 
levels . 
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Figure  3.  Multiple  Shaker  Equivalent  Coincident  (-  -  -)  and 

Quadrature  ( - )  Frequency  Response  Functions  for 

Degree  of  Freedom  2 


We  can  improve  the  mode  isolation  even  further  by  calculating  a  new  set  of 
more  refined  force  levels.  First,  we  form  a  new  truncated  mode  matrix, 


167.539 

165.647 

153.345 


-219.080 
-  23.456 
273.051 


170.043 

-170.031 

109.100 


These  refined  shaker  location  mode  shape  values  were  obtained  from  the  multiple 
excitation  frequency  response  functions  established  with  [F]^.  The  shaker 
location  mode  shape  values  for  the  first  mode  were  established  from  the 
frequency  response  functions  obtained  with  the  first  column  of  [F]j^.  The 
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second  and  third  mode  shape  values  were  obtained  from  the  transfer  functions 
corresponding  to  the  second  and  third  columns  of  [F]^,  respectively.  By 
substituting  the  truncated  mode  matrix  [4»]2  into  Equation  (19)  we  obtain  a 
refined  set  of  force  levels: 


0.624 

-  0.779 

0.615 

1.0 

-  0.138 

-  1.0 

0.586 

1.0 

0.408 

To  obtain  the  frequency  response  functions  associated  with  [Fj^j  we 
repeat  the  scaling  and  summing  operations  as  before.  Note  that  these  opera¬ 
tions  are  to  be  performed  on  the  original,  single-shaker  data  presented  in 
Figure  2.  The  resulting  frequency  response  functions  for  coordinate  are 

also  presented  in  Figure  3  and  are  labeled  (e),  (f),  and  (g).  A  review  of 
these  three  transfer  functions  indicates  that  the  isolation  of  each  target 
mode  is  adequate  and  that  all  the  mode  vectors  can  now  be  established. 

The  determination  of  the  force  levels  needed  for  mode  isolation  only 
involved  the  shaker  location  frequency  response  functions.  Once  the  proper 
force  levels  are  established,  the  frequency  response  functions  associated  with 
coordinates  other  than  the  shaker  location  coordinates  can  be  scaled  and 
summed  as  discussed  above.  This  operation  only  needs  to  be  performed  once 
since  the  iterative  part  of  the  procedure  involves  only  the  shaker  location 
coordinates.  The  mode  shapes  can  now  be  established  by  extracting  from  the 
multi-shaker  equivalent  frequency  response  functions  the  peak  quadrature 
response  at  the  natural  frequency  of  the  mode. 

Having  established  the  modes  in  the  first  isolation  group  we  can  now 

proceed  to  the  second  group.  There  are  two  modes  in  this  group;  therefore,  we 

will  need  two  excitation  locations.  We  select  coordinates  X-  and  X  since 

2.  6 

frequency  response  data  already  exists  for  these  locations.  From  the  transfer 
functions  obtained  with  single-shaker  excitation  at  coordinate  X^  (first 
column  of  Fig.  2)  we  obtain  estimates  of  the  mode  shape  values  for  the  first 
of  the  two  modes  in  the  isolation  group.  For  the  second  mode  in  the  isolation 
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group,  we  obtain  estimates  from  the  frequency  response  functions  associated 

with  single-shaker  excitation  at  coordinate  X,  (third  column  of  Fig.  2). 

o 

Next,  we  form  for  this  isolation  group  the  truncated  mode  matrix 

i.e.  , 


8.585  0.086 

0.137 


and  use  Equation  (19)  to  establish  the  required  force  levels  [F]^.  These 
calculated  forces,  scaled  to  a  maximum  value  of  unity,  are: 


[F] 


1 


1.0  -  0.016 
0.013  1.0 


The  calculated  force  patterns,  for  all  practical  purposes,  are  equivalent 
to  single  shaker  excitation,  and  this  is  how  in  practice  the  modes  should  be 
excited.  However,  for  the  purposes  of  this  simulation  we  will  repeat  the 
frequency  response  function  scaling  and  summing  operations  to  establish  the 
multi-shaker  equivalent  transfer  functions.  From  these  we  can  then  extract 
the  peak  quadrature  response  values  and  establish  the  mode  shapes. 

The  multi-shaker  equivalent  frequency  response  functions  for  coordinates 
and  Xg  are  presented  in  Figure  k.  Initially,  it  would  appear  that  the 
first  three  modes  of  the  system  are  contaminating  the  fourth  and  fifth  modes. 
This  is  correct  if  one  considers  the  coincident  component  of  response  as  the 
measure  of  mode  isolation.  However,  here  we  are  only  interested  in  the 
quadrature  component  when  determining  mode  shapes  and  natural  frequencies. 
Damping  estimates  only  require,  in  addition  to  the  natural  frequency,  the 
frequencies  associated  with  the  peaks  of  the  coincident  response.  As  can  be 
observed  in  Figures  4(a)  and  4(d),  the  frequencies  at  which  these  peaks  occur 
are  not  affected  significantly  by  contamination  from  the  other  modes. 
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Figure  k.  Multiple  Shaker  Equivalent  Coincident  (-  -  -)  and 

Quadrature  ( - )  Frequency  Response  Functions  for 

Degrees  of  Freedom  1  and  8 

The  five  identified  mode  shapes,  natural  frequencies  and  damping  values 
are  compared  to  the  true  values  in  Table  2.  As  can  be  seen,  close  agreement 
exists  between  the  identified  values  and  the  exact  data.  The  mass  weighted 
orthogonality  of  the  five  modes  is  presented  in  Table  3.  The  orthogonality  of 
modes  established  with  single-shaker  excitation  is  presented  in  the  column 
titled  "Initial  Excitation."  The  values  in  the  column  titled  "First  Iteration" 
correspond  to  modes  established  in  the  first  iteration.  The  values  in  the  last 
column  correspond  to  the  first  three  modes  identified  in  the  second  iteration, 
whereas  modes  four  and  five  are  still  the  ones  obtained  in  the  first  iteration. 
As  these  data  indicate,  the  test  success  criterion  of  the  mass  weighted  ortho¬ 
gonality  between  pairs  of  modes  being  less  than  0.10  was  satisfied  after  the 
first  iteration. 
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Table  2.  Comparison  of  Identified  and  Exact  Mode  Shapes 
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VI.  SUMMARY 


An  efficient  direct  measurement  mode  survey  test  procedure  has  been  intro¬ 
duced.  The  procedure  derives  from  the  recognition  that  frequency  response 
functions  obtained  with  single  shakers  can  be  scaled  and  summed  to  yield 
frequency  response  functions  corresponding  to  multiple-shaker  excitation. 

Thus,  the  operations  to  establish  force  levels  needed  to  isolate  each  mode  for 
measurement  can  De  performed  numerically  on  a  small  laboratory  computer/data 
acquisition  system.  The  derived  forces  can  then  be  applied  to  the  structure 
and  the  entire  mode  vector  measured.  Thus,  multiple-shaker  excitation  of  the 
structure  is  not  required  until  shaker  locations  and  force  levels  have  been 
established  numerically. 

The  procedure  introduced  herein  can  also  be  used  to  numerically  identify 
mode  parameters,  from  single-shaker  frequency  response  functions,  with  tradi¬ 
tional  multi-shaker  sine-dwell  testing  techniques.  All  that  is  required  is 
that  frequency  response  functions,  corresponding  to  each  excitation  location, 
be  measured  for  all  accelerometer  locations.  This  quantity  of  data  is  already 
collected  for  many  of  the  existing  analytical  mode  identification  procedures. 
Appropriate  force  levels  can  be  established  to  isolate  each  mode,  and  all 
corresponding  transfer  functions  can  then  be  scaled  and  summed  numerically. 

The  mode  parameters  are  then  established  directly  from  the  resulting  multi¬ 
shaker  equivalent  frequency  response  functions.  It  should  be  noted,  however, 
that  this  approach  still  needs  to  be  tested  with  experimental  data.  Also,  as 
with  any  other  mode  survey  test  procedure,  the  test  article  should  not  be 
released  until  modes  of  acceptable  quality  have  been  established. 

The  numerical  mode  identification  procedure  presented  herein  has,  in  the 
opinion  of  the  author,  several  advantages  over  other  numerical  curve  fitting 
procedures.  However,  the  principal  advantage  is  that  the  dynamic  properties 
are  established  directly  from  scaled,  linear  combinations  of  single-shaker 
frequency  response  functions.  No  analytical  curve  fitting  is  required.  Thus, 
the  disadvantage  of  interpreting  numerical  curve  fitting  results,  as  required 
by  the  existing  post-test  mode  identification  procedures,  is  avoided.  The 
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proposed  procedure  will  not  yield  computational  modes  since  any  mode  to  be 
identified  has  been  selected  by  the  engineer  from  a  review  of  the  raw 
frequency  response  data. 
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